Quantum-classical transition for an analog of double-slit experiment in complex 
collisions: Dynamical decoherence in quantum many-body systems 
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We study coherent superpositions of clockwise and anti-clockwise rotating intermediate complexes 
with overlapping resonances formed in bimolecular chemical reactions. Disintegration of such com- 
plexes represents an analog of famous double-slit experiment. The time for disappearance of the 
interference fringes is estimated from heuristic arguments related to fingerprints of chaotic dynamics 
of a classical counterpart of the coherently rotating complex. Validity of this estimate is confirmed 
numerically for the H+D2 chemical reaction. Thus we demonstrate the quantum-classical transition 
in temporal behavior of highly excited quantum many-body systems in the absence of external noise 
and coupling to an environment. 

PACS numbers: 03.75.-b, 34.10.+X, 05.45.Mt 



The famous double-slit experiment [IJ provides the 
most vivid demonstration of quantum coherent superpo- 
sitions. These manifest themselves in interference fringes 
of the intensity due to the interference between the mat- 
ter waves emerging from different slits. The double-slit 
experiments have successfully demonstrated the wave na- 
ture of, e.g. electrons, neutrons, atoms, small molecules, 
noble gas clusters and fullerenes A process of con- 
verting coherent superpositions into classical sum of in- 
tensities manifests itself in the disappearance of the in- 
terference fringes. This fundamental physical process of 
the emergence of classical dynamics from the quantum- 
mechanical description is referred to as the quantum- 
classical transition (QCT). 

There are two possible roots to describe QCT. The 
first one is decoherence [2| due to the external noise or 
coupling to the environment. This mechanism of QCT 
results from a non-unitary evolution of the system. The 
QCT due to the decoherence for the double-slit experi- 
ment with fullerenes has been demonstrated Q. On the 
contrary, dynamical decoherence [H describes the QCT 
without coupling to environment and only due to the 
intrinsic unitary evolution of a pure quantum state 
This process, unlike decoherence 0, describes the QCT 
on a finite time scale shorter than Heisenberg time, which 
diverges in the macroscopic limit Q • 

In this communication we address the problem of 
quantum-classical transition in the absence of any cou- 
pling to the environment, revealing an effect analogous 
to dynamical decoherence [4|, [5|. We focus on the QCT 
in a temporary quantum evolution. This problem is mo- 
tivated by the work [5j , where time- integration appeared 
to be a precondition for the quantum-classical transi- 
tion. However, without disappearance of the interfer- 
ence fringes at fixed moment of time the QCT is clearly 



incomplete. This rises the question whether dynamical 
decoherence [1 can lead to a QCT in a way that deco- 
herence does 0]. 

Instead of a single-particle problem we consider 
coherent superpositions of clockwise and anti-clockwise 
rotating many-body intermediate complexes (IC) with 
strongly overlapping resonances. Such IC can be cre- 
ated in atomic cluster collisions, bimolecular chemical 
reactions and heavy-ion collisions. The physical pic- 
ture of rotational wave packets and their interference, 
originally revealed for heavy- ion collisions @, 0, HI , has 
been strongly supported by numerical calculations for, 
e.g. the H+D 2 @, F+HD [10] and Hc+H+ llj state- 
to-state chemical reactions. The many-body aspect is of 
primary importance for the QCT since the macroscopic 
world consists of complex systems. We demonstrate here 
that isolated quantum systems can undergo a QCT and 
we obtain the characteristic time for such transition. 

We consider the spin-less collision partners in the en- 
trance (a) and exit (6) channels. The energy fluctuating 
^-matrix elements are taken in the pole form This 
is a good approximation in the regime of overlapping res- 
onances of the IC, Y/D ^> 1, where Y is total width of 
resonance levels and D is average level spacing of the 
IC, provided n c h ^> Y/D with n c h being a number of 
open channels [1^ . The high excitation energy of IC, 
the smallness of the average level spacing of the IC, the 
strong overlap of the resonance levels and the condition 
rich ^ 1 imply that we are deeply in the semi-classical 
region. In this regime we expect the quantum evolution 
to reveal fingerprints of classical dynamics of the classi- 
cal counterpart of the system. We shall see that this is 
indeed the case. 

Under the above stated conditions the time power 
spectrum for the time-delayed collision has the form 
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II(t,0) = H{t)ey^{-t/t u )\V{t,6)\ 2 , where 

w) = ^E( 2J+1 ) M/ ( J ) 1/2ex p( i<i>J )^( ) 



^c^exp(- l ^ 7 t/?i). 



(1) 



tation of the IC uA 



Here, t is the time, is the scattering angle, J is the total 
spin of the IC, and H (t) is the Heaviside step function. 
The deflection angle $ is given by the first J derivative 
of a sum of the potential phase shifts in the entrance 
and exit channels, and the Pj(0) are Legendre polyno- 
mials. In the H(t, 9), the factor exp(—t/tit) accounts 
for the finite life-time, tu = h/T, of the IC. Therefore 
its energy levels have a finite width, h/tu, correspond- 
ing to the continuous spectrum. On the other hand, 
even though the resonance levels are strongly overlap- 
ping, P(t,8) — \V(t,6)\ 2 describes the time evolution in 
discrete spectrum. Note that while P(t, 9) is the prob- 
ability, except for the omitted factor exp(— t/tit), for 
the IC to decay into an angle 9 at the moment i, it 
also describes the time evolution of the angular orien- 
InEq. ©, ci = c M 7 /[RF] 1/2 , 

I T r\\ Ja(b) , J'aCb) 

a - <a / fj, i fx 1 1* ( c m = °)' > and 7^ are 
real partial width amplitudes, and E?, El are the res- 
onance energies. The overlines stand for averaging over 
resonance states. The number of resonance states in- 
cluded in the sum in Eq. CQ) is N = X/D > 1, where the 
energy interval X should be sufficiently long to resolve the 
interference fringes Q. In Eq. Q, W(J) = (|<5S* J (£;)| 2 } 
is the average partial reaction probability taken in the 
Gaussian form, exp[— (J — I) 2 /d 2 ], where / is the aver- 
age total spin of the IC and d is the J-window width. 

The intensity \V(t, 9)\ 2 includes the sum 
jfj^^^i^i ex P[*(^v _ ^D*/^]- We calculate these 
sums by performing first a partial (/x, ^-summation with 
(E'l — E'l ) fixed within an uncertainty of about a few 
units of D. Then, we get instead of the product c^c'l its 
average over resonances, cr[c^' . The latter is given [|| by 
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(l/n)Dp\J- J'\ 



[El - E£ - hcj(J - J')} 2 +(3 2 (J- J') 2 



(2) 



where (3 3> D is the phase relaxation width and ui is the 
angular velocity of the coherent rotation of the highly- 
excited IC. The spin diagonal contribution is evaluated 
by performing first a partial (/i, /i')-summation with fixed 

(p — (j,') value. This yields c^c^, = 5^1. The above ap- 
proximation is justified if the time for formation of the 
many-body eigenstates of the effective Hamiltonian of the 
IC originated from the intramolecular energy relaxation 
due to, e.g. a strong anharmonic coupling between the 
normal vibrational modes, is faster than a characteris- 
tic rotation time of the IC. The diagonal approximation, 
for the states with the same J-values, is the standard 
assumption in the theory of quantum chaotic scatter- 



ing |l5|. The employed approximations imply that in- 
stead of the products c^c^ we can use their average over 
resonances. This substitution does not violate the time 
quasiperiodicity of \V(t, 9)\ 2 , which is characteristic of 
a unitary evolution, since it is still given by a discrete 
superposition of a finite number of Fourier components 
with discrete frequencies (E^ — E'l )/fi. 



We first calculate the sums 
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-£c£cJ'exp[i(£tf -Efrt/h] 

[lis 



(3) 



in the macroscopic limit h — ► 0. Since D oc %* , where 
fO is a number of degrees of freedom, the macroscopic 
limit corresponds to a continuous spectrum approxima- 
tion. Changing from the summation to integration in 
Eq. ([3]) we obtain 



exp[-iut(J - J')} exp[-/3t\ J - J'\/h] 



(4) 



Therefore the continuous spectrum approximation results 
in an irreversible decay of the spin off-diagonal correla- 
tions implying a non-unitary non-periodic time evolution. 
In order to preserve a unitary time-periodic evolution we 
must evaluate the sums ([3]) for a discrete spectrum. We 
do so using first the equidistant spectrum approximation. 
Such an approximation is justified for t <C 2nh/D, when 
the spectrum is not resolved and thus the detailed spec- 
tral properties do not affect the time behavior. Then, for 
X ^> fou),(3, the original sums ((3]) are reduced to 



(l/ir)D0\J-j'\ 



exp(—ikDt/h) 



h— — oo 



with k being integers. We calculate these sums using the 
Poisson summation formula and obtain 

oo 

exp[-27r(j9/£>)|J- J'\\M - Dt/2nh\] 

M=— oo 

exp[i2n(hu;/D)(J - J')(M - Dt/2nh)], (5) 

where M are integers. The above expression is periodic 
in time, with the period 2irh/ D, resulting in a periodic- 
ity of \V(t, 9)\ 2 with this same period. However, we are 
interested in the time evolution of the IC with strongly 
overlapping resonances corresponding to relatively short 
time intervals < t <C 2irh/D. Then, the above sums 
are dominated by the M = term since, for (3 D, 
the M/0 terms are exponentially small and can be ne- 
glected. Then, the M = term equals the expression 
((4]) and we observe that, for t < irh/D, the continuous 
and equidistant spectrum approximation yield identical 
results. Yet, while the continuous spectrum approxima- 
tion violates a unitary periodic in time evolution, the 
equidistant spectrum approximation obviously leads to a 
unitary time-periodic evolution. We have also calculated 
the sums §S§ numerically with E^ and E^ discrete spec- 
tra being uncorrelated for J ^ J' and each having the 
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universal spectral fluctuations of random matrix theory 
(Wigner-Dyson statistics [l6|). We have found that, for 
D <C (3 and t <C 2nh/D, these sums are accurately given 
by the expression (fj|, i.e., by the M — term in Eq. ([5]). 
Again, this is because for t <C 2irh/D the spectrum is not 
resolved and the sums ([3|) are insensitive to the detailed 
spectral properties. Altogether, for t <C 2-kTl/D, we have 

P(t, 9) oc ^(2J + 1)(2J' + 1)[VF(J)VF(J')] 1/2 
Pj{0)P.j'{0) exp[i($ - wt)(J - J') - /?| J - J'\t/h). (6) 

Obviously, the above result can also be obtained within 
the continuous spectrum approximation which is valid 
only for t < Ti/D when the spectrum is not resolved. 
Clearly the appearance of the decaying factor exp(— (3\ J— 
J'\t/h) for t < Ti/D does not violate the unitary charac- 
ter of evolution of the IC. The problem of experimental 
determination of P{t, 9) has been discussed in 0, Il7|. 

We employ the near-side far-side decomposition [11| 
of the collision amplitude, P{t,0) = V^(t,e) + 
V {+ \t,9), so that P(r.,0) oc \V^ +) (t, 0)| 2 + \V^\t, 0)| 2 + 
2Re{P^(t,d)V < -- ) {t,9)*}. The amplitudes V {±) {t,9) 
have the form of Eq. (P) but with Q^ ±} (9) = 
(l/2)[Pj(0) T (2i/ir)Qj(9)} instead of Pj{9). Here Qj{9) 
are Legendre functions of the second kind of degree J and 
the traveling Legendre functions have the semi-classical 
asymptotics Q i j z} (9) ~ [l/^TrJsinfl] 1 / 2 exp{±i[(J + 
l/2)0-7r/4]}. 

It is expected that the quantum superpositions orig- 
inate from the interference between the V^~'(t,6) and 
V^ + \t, 9) in a close vicinity of = 0, n, where the wave 
packets overlap. However, for these angle ranges, the 
asymptotic form of Qj (9) is not valid. Yet, t q ~ c i may 
be conjectured from the time dependence of the near-side 
and far-side intensities, \V^(t, 0)| 2 , for the intermediate 

angle interval, where the asymptotic form of Qy'(9) is 
an accurate approximation. 

We calculate \V {± \t,9)\ 2 using the same approxima- 
tions as those employed for the calculation of P(t,9). 
The \V {± \t,9)\ 2 have the form of Eq. © but with the 

asymtotics of Q^(0) instead of Pj(9). Calculations 
analogous to those of || yield that, for t <C h/d(3, 
\P^(t,9)\ 2 describe Gaussian rotational wave packets 
having a constant width (angular dispersion) ~ 1/d. 
Thus, for t <C h/d(], there is no spreading of the 
wave packets and the memory about their width at 
t = remains. This is characteristic of purely quan- 
tum integrable-like evolution while the classical counter- 
part of the many-body IC is clearly chaotic. On the 
contrary, for t 3> h/d(3, the width of the wave pack- 
ets becomes much greater than that at t = and the 
memory about the angular dispersion at t = is lost. 
In addition, for t > Ti/d(3 and f3t/h > 1, \V i±] (t, 0)] 2 
approach a uniform distribution exponentially fast [6|, 
|:P (±) (M)| 2 ocexp[2cos($±0-wt)/exp(/3t/n)]. We con- 
jecture that this long-time behavior may be considered as 
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FIG. 1: (Color online) Time power spectra related to the 
H+D 2 (ui = 0,ji = 0) -»HD(iij = 3,j f = 0)+D chemical 
reaction for different moments of time. Solid lines correspond 
to d = 2 and (3 = 0.003 eV, dashed lines d = 5 and (3 = 0.003 
eV, and dashed dotted lines d — 10 and (3 = 0.003 eV. 



fingerprint of chaotic dynamics of the classical counter- 
part of the IC. Since the crossover from the integrablc-likc 
to the chaotic-like dynamics of the quantum evolution oc- 
curs between t\ <^Ti/(3d and S> h/j3d, the characteris- 
tic time for the QCT can be evaluated as £ g _ c ; ~ h/ftd. 
But then another feature of the QCT should be the dis- 
appearance of the interference between waves originating 
from the decay of rotating in opposite directions IC. This 
transformation of the quantum superposition into an in- 
coherent sum of the intensities should also happen at 
tq-d ~ h//3d. 

To test this we use Eq. ^ to describe rotational wave 
packets and their interference revealed [9( for the chem- 
ical reaction H+D 2 (wi = 0,ji = 0) — >HD(u/ = 3, j/ = 
0)+D, where Vi,ji and Vf, jf are vibrational and rota- 
tional quantum numbers for the initial and final states, 
respectively. In this reaction, the wave packets likely 
originate from interference of the overlapping resonances 
of the IC fl8| . From the numerical calculations Q we 
deduce <& = 135° and hid = 0.045 eV. It was found that 
the main contribution to the time-delayed reaction mech- 
anism comes from the total angular momenta J = 15—20. 
Accordingly, we choose I = 18. To reproduce the width 
of the wave packets, ~ 25°, at an early stage, when ~ $, 
we take d = 2. To account for the strong interference 
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FIG. 2: (Color online) Classical limits for the time power 
spectra presented in Fig. 1 (see text). 



contrast in forward direction we choose (3 — 0.003 eV. 

In Fig. ED we present P(t,9)/P dlag (9), where P dla g(9) 
is a time-independent quantity given by the diagonal 
J = J 1 terms in Eq. ([6]). The reason for this scaling is 
that amplitudes of individual wave packets are strongly 
and abruptly enhanced in the vicinity of 9 = 0, 7r due 
to the azimuthal symmetry of the problem [8J. How- 
ever, Pdiag(0) has the similar enhancement for 9 ~ 0, n, 
so that the quantity P(t, 9) / Pdiag{9) permits us to probe 
interference between the wave packets whose amplitudes 
depend smoothly on 9. Another advantage of this scal- 
ing is that Pdiag (9) corresponds to the very quick phase- 
relaxation (Ti/(3 — ► 0) yielding the universal limit of 
the theory of quantum chaotic scattering In this 

limit, wave packets are uniformly spread and their in- 
terference is completely destroyed due to the absence of 
the spin off-diagonal contributions. Thus, a deviation of 
P(t,9)/P ( iiag{9) from a constant unity is a quantitative 
measure of the deviation of the collision process from the 
universal limit of the quantum chaotic scattering theory. 

In Fig. [1] we present P(t, 9)/ Pdiag{9) for three moments 
of time. One can see that the wave packets, initially ori- 
ented along 9 — <£>, 27r — $, rotate towards each other. 
Interference fringes are produced when the wave pack- 
ets start to overlap. For d = 2, (3 = 0.003 eV we have 
t q - c i ~ Ti/ (3d — 220 fs. This time is longer by about 
a factor of 6.5 than the longest time 3T/8 in Fig. [TJ 
where T — 2tt/lli — 90 fs. We observe a strong con- 



trast of the interference fringes. Now we decrease t q - c i 
by taking d — 5 and keeping (3 the same. Then, for 
t = 3T/8 ~ 0.38 x t q - c i, the interference fringes are sup- 
pressed but still visible in the bottom panel of Fig.[TJ We 
further decrease t q _ c i by taking d — 10 with the value of 
(3 unchanged. Then, for t = 3T/8 ~ 0.77 x t g _ d , the 
interference fringes practically disappear manifesting a 
transformation of the coherent superpositions into the 
incoherent mixture of the intensities. Thus, we have con- 
firmed the validity of our estimate for the t q — c i. This 
estimate has also been confirmed by calculations where 
we have increased (3 keeping d = 2 fixed. 

To test that the interference fringes do originate from 
the interference between the P M (t,9) and V ( ^(t 7 9) 
we calculate p(+)(i,0) + P^(t,9), where P^(t,9) = 
\V^{t, 9)\ 2 . We employ the same approximations used 
in calculating P(t,9), Eq. ©. We obtain P^(t,9) hav- 
ing the same form as Eq. © but with (9)Q ( ^\9)* 
instead of Pj(9)Pj,(9). In Fig. EU we present [P(+)(t, 0) + 

P^(t,9)}/[P^J g (t,9) + P$ g (t,0)] for the same values 
of parameters as used in Fig. ED The time independent 
quantity [Pj^ g (t,9) + PjfagfaO)] contains the diagonal 
J = J' terms only. The reason for this scaling is the same 
as that employed in Fig. EQ In addition, Q J ± \9) have a 
logarithmic divergence at 9 = 0, tt resulting in unphysi- 
cal increase of (t, 9) in a close vicinity of 9 = 0, tt. 
The presentation in Fig. fallows to scale out these unde- 
sirable features enabling a meaningful comparison with 
Fig. [Q We see that the interference fringes in Fig. [JJ 
do originate from the interference between the V* + '(t, 9) 
and V^~\t,9) while Fig. [5] shows a classical sum of the 
near-side and far-side intensities. 

Note that the linear size of the chaotic billiard is 
only about twice as bigger as the spatial width of the 
wave packet at t = 0. Therefore, at any moment of time, 
the width of the wave packet can not exceed its width at 
t = more than by about a factor 2. From the obtained 
criterion the interference fringes wash out when the width 
of the spreading wave packets begins to exceed its initial 
value at t — by a factor > 2. It could be of interest to 
check if the interference fringes, at fixed moments of time, 
for chaotic billiard would disappear upon increasing the 
size of the billiard by a factor of 2-3 keeping the width of 
the wave packet at t = the same as in [f|. 

Our approach does not intend to substitute rigorous 
numerical calculations 0]. Yet, it does show that the 
complicated many-body collision problem can be accu- 
rately represented by a much simpler physical picture of 
a weakly damped (/3 <C TvjS) quantum rotator. This map- 
ping has enabled us to address a problem of the quantum- 
classical transition in complex collisions and evaluate a 
characteristic time for such a transition. Remarkably, 
this simple picture of a weakly damped coherent rotation 
has also emerged from the numerical calculations for the 
F+HD [l(| and He+Hj [ll[ chemical reactions. Origi- 
nally this same physical picture of stable rotational wave 
packets was revealed for heavy- ion collisions @, Qj| ■ In- 
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terestingly, this is in spite of the fact that a characteristic 
rotation time of the IC formed in bimolccular reactions 
is ^8 orders of magnitude longer than that for the IC 
formed in heavy-ion collisions. 

The obtained criterion yields i 9 _ c / ~ t spr /d, where 
t S pr = fi//3 is a characteristic time for the spread of a 
wave packet [6]. For classical macroscopic systems the 
wave packet spreading may be associated with the chaotic 
dynamics and t spr could be macroscopically large. How- 
ever, since d is in h units and, in the macroscopic limit, 
h vanishes, then t q _ cl oc % — > 0. This provides an alter- 



native to the explanation in [2| of the practically instan- 
taneous collapse of coherent superpositions into classical 
sum of intensities for macroscopic systems. Unlike deco- 
herence , in our approach "the quantum origins of the 
classical" are obtained in the absence of external noise or 
coupling to an environment and only due to the unitary 
evolution of a pure quantum state. 
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